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In order to stabilize Majorana excitations within vortices of proximity induced topological superconductors, 
it is mandatory that the Dirac point matches the Fermi level rather exactly, such that the conventionally 
confined states within the vortex are well separated from the Majorana-type excitation. Here, we show 
by angle resolved photoelectron spectroscopy that (Bii_a^Sba,)2Te3 thin films with x = 0.94 prepared by 
molecular beam epitaxy and transferred in ultrahigh vacuum from the molecular beam epitaxy system to the 
photoemission setup matches this condition. The Dirac point is within 10 meV around the Fermi level and 
we do not observe any bulk bands intersecting the Fermi level. 

A topological insulator (TI) is characterized by a 
bulk energy gap which hosts conducting helical sur¬ 
face states.These surface states are protected by 
time reversal symmetry according to a topological 
Z2 invariant.^’^ A non-trivial Z2 number, implying 
helical surface states, is induced, e.g., by spin-orbit 
coupling which inverts the band order around the gap 
at some high symmetry points of the Brillouin zone.^~^ 

The resulting non-degenerate surface bands are the 
starting point for the search of Majorana excitations 
(ME’s).^ These excitations being its own antiparticle 
are rigidly pinned to the Fermi level They are 

discussed, e.g., as a possible pathway towards topolog¬ 
ical quantum computation^^’using its non-Abelian 
braiding statistics and its intrinsic high degeneracy. 

Several tunneling experiments have found signatures 
of MB’s in solid state systems, such as at the end of 
a semiconducting nanowire in contact with a super¬ 
conductor within an axial magnetic field,at the 
end of ferromagnetic chains on top of a superconductor 
with spin-orbit coupling,at the surface of possible 
topological superconductors^^”^^ or at the surface of 
superconductors on top of topological insulators. 

However, none of these experiments are sufficiently 
controlled to exclude all other possible explanations^^”^^ 
for the found peaks at Bp. Thus, a smoking-gun 
experiment is still mandatory. 

Such an experiment could use a TI with a pierced su¬ 
perconductor on top probed by scanning tunneling mi¬ 
croscopy (STM).^^ The pierced holes will act as traps 
for the quantized vortices induced by magnetic field B. 
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If the number of vortices 
within a hole is odd, it will 
provide a ME at 
with predicted spatial 
distribution,^^ while, if 
even, at slightly higher or 
slightly lower B, the ME 
will disappear in favor of 
conventional excitations away 
from Ep^. In order to realize 
such an experiment, the MB’s 
have to be energetically well 
separated from conventional 
confined states within the 
topological vortex, i.e. the 
density of states of the Dirac 
cone of the TI has to be 
sufficiently low.^^ If the hole 
size is about the coherence 
length of the superconductor, one can estimate the 
energetic distance of the ME to the f irst excited state by 
El — Ep A^/+ (Ep — Ed)^ with Ep being the 
Dirac point energy and A the superconducting gap.^^’^^ 
Eor a typical energy resolution of O.I meV at T = 0.3 
and A ^ I — 2 meV, this implies |Ed — E^p| < 10 meV. 
Additionally, the Dirac point has to be energetically 
remote from any bulk states. 

This cannot be achieved with binary TFs as, e.g., 
Bi2Te3, Sb2Te3 or Bi2Se3.^^ Eor Bi2Te3, Ep is buried in 
the bulk valence band (BVB) and Ep is located in the 
bulk conduction band (BCB), whereas for Sb2Te3, Ep 
is in the bulk energy gap and Ep is located in the BVB.^^ 

Mixing these two compounds leads to the ternary 
system (Bii_a,Sba,)2Te3, which exhibits the same 
tetradymite structure as Bi2Te3 and Sb2Te3 with a mix 



FIG. 1. Atomic struc¬ 
ture of (Bii_a,Sba;)2Te3 
(tetradymite-type crys¬ 
tal) with colors of differ¬ 
ent atoms marked. 
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FIG. 2. (a) XRD 20/0 scan; (000/) reflections of the grown 
thin film and Si(lll) substrate reflections are labeled; (b) 
RBS measurement to determine the stoichiometry; Sb Te 
peak and Bi peak are labeled. 


of Sb and Bi atoms in one of the layers (Fig. 1 ). Previous 
angle resolved photoelectron spectroscopy (ARPES) 
and transport experiments have shown that this ternary 
alloy exhibits Ey — E’d — 70 meV at x ^ 0.94 and 
Ey — Ey> ^ — 30 meV at x « 0.96.^^ This indicates that 
the favorable condition of Ey — Ey> is possible at an x in 
between. Moreover, the transport data are compatible 
with a charge carrier concentration originating from 
the Dirac cone only, such that one can anticipate the 
absence of bulk bands at Ed, albeit the authors did not 
study the bulk bands in detail by ARPES. 

Here, we investigate thin films of (Bii_a,Sba,) 2 Te 3 on 
Si(lll) by ARPES, after an ultra-high vacuum (UHV) 
transfer from the molecular beam epitaxy (MBE) cham¬ 
ber to the ARPES setup. We established a sufficient 
accuracy of |Ed - Ef| for the selection of samples for 
the mentioned Major ana experiments, which is applied 
to several stoichiometries {x = 0.48, 0.82, 0.93 and 
0.94). Eor X = 0.94, we show that |Ed - Ef| < 2 meV 
±7 meV, which has to be contrasted with the previous 
precision of about 50 meV of [38,40]. Eurthermore, 
we confirm the absence of bulk states at Ep, at least, 
for photon energies hv = 8.44 eV and 21.2 eV. A 
similar condition has been achieved previously for the 
quaternary compound BiSbTeSe 2 ,^^’^^ but not for any 
ternary compound. Such a ternary might be more 
easy to handle in terms of stoichiometry and interface 
chemistry, which is probably also favorable for transport 
and spin transport experiments requiring Ed close to 
Ep.38,41-43 ^0 show the data for x = 0.94, 

where Ed is closest to Ep. 

Thin films of (Bii_a,Sba,) 2 Te 3 were grown by MBE 
on highly doped Si(lll) substrates, prepared in-situ by 
flash annealing in UHV. Knudsen effusion cells were used 
for the evaporation at a substrate temperature of T = 
275°C. X-ray diffraction (XRD) and reflectivity (XRR) 
measurements were performed ex-situ. The consecutive 
order of reflections correspond to ( 000 /) reflections of 
hexagonal Bi 2 Te 3 , evidencing that the (Bi,Sb) 2 Te 3 film 
is of single crystal nature (Eig. 2(a)). XRR reveals 
a film thickness of 23 ± 1 nm. Rutherford backscat- 
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FIG. 3. ARPES data of (Bio.o 6 Sbo. 94 ) 2 Te 3 , hn — 21.2 eV, T— 
200 K (MBE growth, in-situ transfer): (a),(b) E(A;||) disper¬ 
sion along r - K at different resolution (black: high intensity, 
white: low intensity), Fermi level Ep as determined on poly¬ 
crystalline Gu is marked as a solid black line; (c) Momentum 
distribution curves (MDGs) between binding energies Eb = 
0.240 eV and Eb = -0.110 eV, continuously offset for clarity; 
red and black MDG’s are at Eb = 0.170 eV and Eb = 0.045 
eV, respectively; green MDG corresponds to Fig. 4(a) 


tering spectrometry (RBS) determies the Bi yield to 
2.10 X 10^^ atoms/cm^ and the merged Sb and Te yield 
to 8.8 X 10^^ atoms/cm^ (Fig. 2(b)). Combined with 
the XRD measurement, the stoichiometry is calculated 
to be X = 0.94 ± 0.01 (for a more detailed discussion, see 
supplements^). 

The samples were transferred in-situ from the MBE 
system to the ARPES system with a vacuum shuttle at 
pressure p = 5 x 10“^^ mbar. This prevents oxidation 
and surface contamination of the sample. The ARPES 
spectra were recorded at T = 50 K with a MBS Al hemi¬ 
spherical analyzer providing an energy resolution of 10 
meV and an angular resolution of 0.1°. Two different hn 
are used provided by the non-monochromatized He la 
resonance {hn = 21.22 eV) and monochromatized light 
from a microwave-driven Xenon source {hn = 8.44 eV). 
Eor both hu, Ey has been determined independently on 
polycrystalline Cu with an accuracy of 7 meV. There are 
no indications of time-dependent band bending within 2 
days in contrast to similar experiments on Bi 2 Se 3 ,^^’^^ i.e. 
the band bending is below 7 meV. Photovoltage effects 
on (Bio.o 6 Sbo. 94 ) 2 Te 3 have been checked to be negligible 
by the temperature dependence of the spectra between 
50 K and 300 K. The directions in momentum (/c||) space 
of the ARPES spectra are deduced by comparison with 
ab-inito calculations for Sb 2 Te 3 ( 0001 ),'^^ i.e., the lips of 
the star-like structure in Eig. 5(a) point to F-M direction. 

The measured band structure in F-K direction is shown 
in Fig. 3(a) —(b). The characteristic linear band disper- 
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FIG. 4. (a) MDC’s (purple dots) at indicated with two 
Voigt fit curves (black, red) with indicated peak positions 
(Max 1, Max 2) and resulting fit (blue curve); (b) same as 
(a) with only one Voigt fit curve (blue) and marked FWHM; 

(c) peak positions deduced from fitting as shown in (a), as 
function of only region 1 is used for fits and the linear 
regression (black, red dashed line) to determine F^d; error 
bars (blue ellipses) result from the fitting procedure; vertical 
error bar (blue) at F^d is the error on the intersection point; 

(d) FWHM’s of the MDC peak close to F^f determined as 
marked in (b). 


sion of the Dirac type topological surface state (TSS) is 
visible (Fig. 3(b)) having the Dirac point close to F^p- 
This can also be deduced from the momentum distri¬ 
bution curves (MDC’s) in Fig. 3(c). To evaluate the 
position of F^d more precisely, two Lorentzian peaks con¬ 
voluted with a Gaussian function (Voigt function^^~^^) 
were fitted to the MDC’s of the TSS (Fig. 4(a)) for 
binding energies F^b = 0.170 — 0.045 eV (region 1, 9 fit 
parameters per MDC, see supplements^). The posi¬ 
tions of the resulting peak maxima are plotted in Fig. 
4(c) exhibiting an average error of Ak = ±1.5 x 10“^ 
A . For F^b = —0.030 — 0.045 eV, the two maxima are 
too close to be separately visible in the MDC’s resulting 
in much larger errors up to Ak = ±1.5 x 10-2 A . To 
determine the position of F^d, a linear fit of the data for 
F^b = 0.170 — 0.045 eV (region 1 of Fig. 4(c)) is extrapo¬ 
lated to F^f- The linear regression, which includes the k\\ 
errors from the Voigt fits and the energy errors resulting 
mostly from the uncertainty of F^p, reveals a Dirac veloc¬ 
ity for the TSS at —k\\ of ud = (3.8 ± 0.2) x 10^ ™ and 
at k\\ of ud = (3.9 ±0.2) x 10^ The intersection of the 
two TSS is at F^b = 2 ± 7 meV. 

Another way to determine F^d is given by the full width 
at half maximum (FWHM) of a single Voigt peak to the 
MDC’s close to F^p (Fig. 4(b)). Due to thermal exci¬ 
tations, this peak can also be recorded above F^p- The 
FWHM’s of the peaks are plotted in Fig. 4(d) showing a 
minimum at F^b = 0 ± 10 meV. Thus, (Bio.o 6 Sbo. 94 ) 2 Te 3 
exhibits F^d = with a precision better than 10 meV. 
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FIG. 5. ARPES data probing the bulk valence band (BVB) 
(MBE growth, in-situ transfer): (a) GEG at Eb = 120 meV 
showing the star shaped maximum area of the BVB, where 
kx and ky are in the F - K and F - M directions, respec¬ 
tively; hu = 8.44 eV; cuts mark the k\\ directions of the plots 
in (b) —(f) as labeled; (b) —(f) ARPES data in E{k\\) repre¬ 
sentation along the cuts marked in (a) with photon energies 
indicated; the solid black line Ef is determined on poly crys¬ 
talline Gu; (g)-(j) GEG’s at (g) Eb =0 meV, (h) 10 meV, (i) 
20 meV, (j) 40 meV, hu = 21.2 eV. 


This is exactly the precondition for the search of ME’s 
by STM as described above, albeit the superconductor 
on top might shift F^d again. 

Another requirement is, that the BVB is further away 
from F^p than the induced superconducting gap. In 
order to probe this, we plot ARPES data through the 
maximum of the BVB in k\\ space (along cut 1 and cut 
2 in Fig. 5(a)) and the TSS (along cut 3 in Fig. 5(a)) 
at two different kz, i.e. two different hu (Fig. 5(b)-(f)). 
The valence band maximum is found at F^b — 60 meV 
for both Fz/, indicating little dispersion along kz as 
expected for the layered material. The absence of the 
BVB at F^f is corroborated by the constant energy 
cuts (CEC’s) at hu = 21.2 eV (Fig. 5(g)-(j)) which ex¬ 
hibit the onset of a BVB band structure at F^b — 40 meV. 
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In summary, we synthesized a ternary topological insu¬ 
lator (Bio.o 6 Sbo. 94 ) 2 Te 3 , which exhibits E'd — = 2 ± 7 

meV and does not show any indications of the bulk 
valence band close to This is favorable for both, 

transport measurements, in particular, probing spin- 
transport^^ and a good starting point for the search of 
Majorana quasiparticles by STM.^ 
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Supplement: Tuning the Dirac point to the Fermi level 
in the ternary topological insulator (Bii_2;Sba;)2Te3 

A. Fitting routine 



FIG. S 1. (a) MDC’s (purple dots) at indicated Eb with two Voigt fit curves (black, red) and 
resulting fit (blue curve); Lorentz and Gauss contributions for the first peak (light green, dark 
green curves) and the second peak (light brown, dark brown curves) are added; (b) peak positions 
deduced from fitting as shown in (a) for Eb — —0.01 eV — 0.24 eV with different regions as 
described in the text marked. 


To quantitatively determine the Dirac point position, the MDC’s of the TSS are fitted, 
firstly, with two Voigt functions. A Voigt function Vi{x) is a convolution of a Lorentzian 
and a Gaussian function, such that: 
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where y ist the offset, Ai the amplitudes, Xc^ the centers of the Voigt functions, /r, the 
profile shape factors and icl, 'Wq the widths of the Lorentzian and the Gaussian contribution, 
respectively. The two Voigt functions were adjusted to the MDC’s (Fig. Sl(a)) within the 
energy range = —0.01 eV — 0.24 eV. In total, the ht has 9 free parameters. The fit is 
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adjusted by a Levenberg-Marquardt algorithm and the standard deviation is read out for 
every parameter. Note that the widths wl and wq are not assigned with an index i, thus, 
the widths are the same for the two peaks at —k\\ and k\\. 

The resulting total amplitude, as well as the Lorentzian and Gaussian amplitudes are higher 
for all —fc||. This is most likely due to the geometrical setup of the ARPES system, i.e., 
matrix element effects induce the difference in intensity. 

The ht results for the peak positions are shown in Fig. Sl(b) including error bars. In the 
region labelled region 2, the error bars are quite large, since the MDC exhibits only one max¬ 
imum. Moreover, the dispersion deviates from linearity for Eb > 0.170 eV. Consequently, 
we used only the energy region Eb = 0.170 eV — 0.045 eV (region 1) for the determination 
of Eb as shown in Fig. 4(c) of the main text. 
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B. RBS and XRD measurement 



Channel 


(b) 




c 10000- 

3 

ifHi 

1000; 

> 

lOOl 

'(/) 

c 

o 

10i 

c 

1, 


0 



io 20 30 40 50 60 70 


26 [degree] 


FIG. S 2. (a) RBS measurement to determine the stoichiometry; Sb Te peak and Bi peak are 
labeled; (b) XRD 2©/0 scan; (OOOZ) reflections of the grown thin film and Si(lll) substrate 
reflections are labeled. 

Rutherford backscattering spectrometry (RBS) using 1.4 MeV He"*" ions at a scatter¬ 
ing angle of 170 degrees was applied in order to determine the stoichiometry of the film. 
Fig. S2(a) shows a typical RBS spectrum. The absolute areal coverage was determined by 
integration of the two peaks. The peak of Bi is clearly separated and yields 2.10 x 10^® 
atoms/cm^. Due to the low difference of the atomic mass of Sb and Te these two peaks 
cannot be resolved separately. Therefore, only the sum of Sb and Te atoms of 8.8 x 10^® 
atoms/cm^ is determined. The experimental accuracy of the absolute values is estimated to 
be i2%. 

X-ray diffraction (XRD) measurements were accomplished employing a Bruker D8 high- 
resolution diffractometer. The XRD scans were carried out in a symmetric 20/0 configura¬ 
tion. Figure S2(b) depicts the XRD curve of the investigated (Bi,Sb) 2 Te 3 sample. Numerous 
peaks originating from the (Bi,Sb) 2 Te 3 are seen, which can be assigned to reflections in the 
(OOOZ) direction with an accuracy of ±0.5°. The consecutive order of reflections correspond 
to (OOOZ) reflections of hexagonal Bi 2 Te 3 , evidencing that the (Bi,Sb) 2 Te 3 film has separate 
BiSb and Te layers according to the stacking sequence shown in Fig. 1 of the main text. 
If the deviation from the crystal stoichiometry (Bi,Sb) 2 Te 3 would be more than ±1%, the 
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resulting distortion of the crystal would shift the XRD peaks by more than ±1°. Thus the 
accuracy of the stoichiometry is estimated to be ±1%. 

Combining the RBS and XRD results, the Sb concentration can be computed to a; = 
94% ± 1% within the formula (Bii_a;Sba;) 2 Te 3 . 
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